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Single site height probabilities in the Abelian sandpile model, and the corresponding mean 
height (h), are directly related to the probability P re t that a loop erased random walk passes 
through a nearest neighbour of the starting site (return probability). The exact values of 
these quantities on the square lattice have been conjectured, in particular (h) =25/8 and 
■Pret = 5/16. We provide a rigourous proof of this conjecture by using a local monomer-dimer 
formulation of these questions. 

I. INTRODUCTION 

As different as they may appear at first sight, the dimer model, the loop-erased random walk 
(LERW) and the Abelian sandpile model (ASM) are very closely connected, as they all have an 
alternative formulation in terms of a unifying concept: spanning trees. 

Among the three models just mentioned, the dimer model is the oldest one as it was formulated 
in 1937 by Fowler and Rushbrooke [l|, although the first exact results were obtained somewhat 
later, in the sixties, by Kasteleyn, Fisher, Temperley and Stephenson j-[6(]. The dimer (or domino) 
model, also known as the perfect matching problem, has been the subject of an increasing number 
of works over the last decade, see for instance [7] and the references therein. The correspondence 
between packed dimer configurations (no vacancy or monomer) on the square lattice and spanning 
trees on a sublattice was established by Temperley [8], and generalized to general planar graphs 
by Burton and Pemantle [9]. This correspondence holds in the presence of monomers, and leads 



id, 



1 if ] . It will be recalled in Section 3. 
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to spanning webs rather than spanning trees 

The loop-erased random walk was defined by Lawler 
symmetric random walk by removing the loops as they appear. 
LERW and the spanning trees has been discussed in several works 
the probability measure on LERW sample paths coincides with the uniform measure on chemical 
paths of spanning trees. 

Finally the Abelian sand pile model is an open stochastic dynamical system, which works like a 
non-linear diffusion process (id. [lij]. The ASM is defined in terms of height variables, attached to 
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the sites of a two-dimensional square lattice. The height variables hi take the four integer values 
1, 2, 3 and 4, but only a fraction of all height configurations, called recurrent, keep reoccuring 
when the dynamics is run over long periods. The long time behaviour of the model is controlled 
by the stationary measure on the height configurations. This measure is uniform on the recurrent 
configurations and vanishes on the non-recurrent ones. The connection to spanning trees stems 
from the burning algorithm [^(J, which establishes a direct and one-to-one correspondence between 
the recurrent configurations of the sandpile on a lattice and spanning trees on the same lattice. 

All three models have non-local features. This is manifest for the dimer model and the LERW, 
whereas the non-locality in the ASM comes from the recurrence criterion [n3], which requires to 
scan the whole of a height configuration before one can declare it recurrent. These non-local 
features usually make explicit calculations particularly hard. Let us also mention that the question 
of conformal invariance has been addressed in these threejnodels, directl y in ter ms of a conformal 
field theory for the dimer model 



spanning webs [24| and the ASM [2514281] , and as stochastic 



Loewner processes (SLE) in the case of LERW and spanning trees 



29|,|30|]. 



The purpose of this work is to compute the value of the return probability for the LERW 



-fret = announced in |31[ |. This value was based on the observation that the return probability 
and the ASM mean height are related by P re t = 2 — I> anc ^ on an earlier conjecture made in 
[27I for the mean height (h) = in the limit of the infinite square lattice. It turns out that 
these two numbers are themselves related to three others, recently introduced by Levine and Peres 
[32] ]. namely the looping constant £, the ratio r between the number of spanning unicycles and the 
number of spanning trees, and the mean length A of the cycle in a spanning unicycle, all defined as 
limits over increasing finite square grids. In addition these five numbers have all a d-dimensional 
analogue, and remain rationally related in any dimension [3^]. In this sense, our result provides a 
proof for the values of these five numbers in two dimensions, P re t = jg, (h) = ^f>£ = f> r = | an d 
A = 8. 

In Section 2, we recall the expressions of the LERW return probability and of the ASM mean 
height in terms of spanning trees with some specific properties. Section 3 contains the proof 
itself and reduces the counting of the required spanning trees to certain local configurations of 
monomers and dimers. We should emphasize that the proof does not rely on the exact evaluation 
of the multiple integral on which the conjecture made in [2^] is based. On the contrary the proof is 
essentially combinatorial, and shows that the counting of spanning graphs with certain non-local 
properties can be reduced to the counting of local monomer-dimer arrangements, which can then 
be easily carried out. It thus avoids the full complexity of the graph theoretical computations, 
inherent in 
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33J . Our proof shows and explains why the LERW return probability and the other 



four related quantities are such simple numbers. 
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II. THE LERW RETURN PROBABILITY AND THE ASM MEAN HEIGHT 



The Abelian sandpile model [lg] in finite volume is defined by height variables, located at the 
sites of a finite, two-dimensional square grid and taking the values 1, 2, 3, 4 in stable configurations. 
Particles are added one by one at a random site, thereby increasing by 1 the value of the height 
at that site. If the height exceeds 4, then the site becomes unstable and topples, transferring 
one particle to each of its neighbouring sites (sand may fall off the system at boundaries). As 
a consequence, one (or more) neighbour may become unstable, in which case it topples too, and 
so on for the neighbours of the neighbours. When no unstable site remains, another particle is 
dropped at a random site and the relaxation process repeated. 

The analysis of this discrete dynamics was performed by Dhar [lj]], who introduced the notion 
of recurrent configuration. He showed that the system enters the subset of recurrent configurations 
after a finite time, which depends on the initial configuration, and never leaves it. Therefore the 
probability distribution which controls the asymptotic behaviour of the sandpile vanishes on non- 
recurrent configurations; it can be shown to be uniform on the recurrent subset (of size ~ 3.21 l2 
for a L x L grid). Dhar also gave a criterion to select the recurrent configurations. Its explicit form 
will not be important for what follows. 

As sandpile configurations are made of random variables valued in {1,2,3,4} and distributed 
according to the uniform measure on the recurrent subset, it is natural to ask about the statistical 
properties of these variables and their spatial correlations. In particular the first question concerns 
the distribution of the height at a single site. As a first step, and in order to avoid boundary effects, 
we consider a site deep in the middle of the grid, and take the infinite volume limit. We denote by 
Pi, for % = 1,2,3,4, the resulting probabilities, namely Pi is the probability that the height at any 
fixed site be equal to i, in the limit of an infinite grid. 

The first of these four probabilities has been obtained by Majumdar and Dhar [34], 
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P\ = — t: =■ ~ 0.07363. (2.1) 

TT Z IT 6 

The other three probabilities P2, P3 and P4 turned out to be more complicated, and are more 
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expressed in terms of spanning trees. As mentioned in the Introduction, the burning 



convenient 

algorithm 20] , which is the algorithmic translation of the recurrence criterion found by Dhar 1_9(] , 
establishes a one-to-one mapping between the set of recurrent configurations and the set of rooted 
(oriented) spanning trees on the same lattice. The characterization of those spanning trees which 
correspond, under this mapping, to recurrent configurations with certain height values at certain 
positions, has been worked out in 33f] . As a result, the fractions Pi of recurrent configurations 



which have a height equal to i at a reference site, are given by the following fractions among 
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spanning trees, 

pi=^, p *= p i+ik> Ps=P2+ ^' P4=P3+ f ' (2 - 2) 

where X^, for = 0, 1, 2, 3, is the number of spanning trees such that the reference site has exactly 



k predecessors among its four nearest neighbours [33J (a site x is called a predecessor of y if the 
path along the tree from x to the root passes through y). The quantity N is the total number of 
spanning trees on the grid. The identities (|2.2|) are valid for a finite grid. Their infinite volume 
limits exist, and define the sought probabilities. 

Unlike Xq, the quantities X/. for k > 1 are more complicated to compute because one has to 
count the spanning trees satisfying a non-local constraint (those in Xq satisfy local constraints). 



These numbers have been first computed in 



331 ] . the results taking the following form, 
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P2 =2-2^ + ^ + 4< ^ 
1 3 1 12 h 3I 2 , tS 
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where Ji and I 2 are two complicated, multiple integrals. Their numerical evaluation gave the 
following values for the probabilities, P2 — 0.1739, P3 ~ 0.3063 and P4 ~ 0.4461. 

Somewhat later, the calculation of Pi on the plane (and the upper half-plane) was reconsidered 
in [2^], which led to an exact relation between Pi and P3, 

3 12 48 

(tt - 8)P 2 + 2(vr - 2)P 3 = tt - 2 - - + — - — , (2.6) 

7T 1T Z 7T 

and a similar relation between the two integrals I\ and l<i. A conjecture on the value of the 
remaining integral, based on its numerical evaluation to twelve decimal places (later pushed to 25 
places), then yielded the following conjectural values for the probabilities [271 ]. 
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and for the stationary mean height (particle density), 



25 

(h) = A + 2P 2 + 3P 3 + 4P 4 = — • (2i 



Conjectured by Grassberger 



351 ] almost 20 years ago, the value 25/8 has remained an enigmatic fact 
of the Abelian sandpile for its striking simplicity. We will prove this conjecture, not by computing 
the mean height within the sandpile model, but rather by using a recent observation concerning 
the LERW return probability. 

For a LERW starting from the origin, we let the return probability P re t be the probability that 
a path passes through a fixed nearest neighbour of the origin, say the right neighbour. In terms 
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of uniformly distributed spanning trees, P re t is the probability that the origin is a predecessor of 
its right neighbouring site. By listing explicitly the situations where this is the case among those 
where the right neighbour has k predecessors among its nearest neighbours, the following identity 



follows 



ail 



" ret = W + W + W (2 - 9) 

By using the equations (|2.2p to trade the for the Pi (and the identity P± + P2 + P3 + P4 = 1), 

the previous combination turns out to be directly related to the mean height of the sandpile model, 

1 15 (h) 5 

P vet = -(-3Pi -P 2 + P 3 + 3P 4 ) = 2 ( p i + 2p 2 + 3^3 + 4P 4 ) - j = Y " 4- ( 2 - 10 ) 

In order to actually compute the return probability, we express directly the fact that the origin 
is a predecessor of its right neighbour, without paying attention to the other nearest neighbours. 
This is pictured in Fig.[H where all possible configurations are shown. 




C C A B 

FIG. 1: Decomposition of th return probability in terms of chemical paths between two nearest neighbours. 

If A, B, C denote the corresponding fractions of spanning trees with the chemical paths as 
shown, we have that 

P ret = ^ + 2(A + B + C). (2.11) 

It remains to compute A, B and C. This will be done in the next section by using particular local 
arrangements of monomers and dimers. 

III. MONOMER-DIMER COMPUTATION OF LERW RETURN PROBABILITY 

Consider a (2n — 1) x (2n — 1) square lattice C with the rightmost lower site removed. The 
sites of the lattice can be subdivided into three subsets: (1) black sites forming the sublattice B 



6 



of sites with odd-odd coordinates; (2) white sites forming the sublattice W of sites with even-even 
coordinates; (3) the other remaining sites colored grey. The corner site removed is black and 
denoted below by r (for root). 

Consider a dense packed dimer configuration on C. Each dimer covers two sites, of colors black 
and grey, or white and grey. Temperley's correspondence associates two sets of arrows to the dimer 
configuration in the following way. We replace each dimer covering a black site by an arrow directed 
from the black site to the grey one, and similarly we replace each dimer covering a white site by 
an arrow directed from the white site to the grey one. The sets of black and white arrows are two 
acyclic configurations of arrows, which form spanning trees on the two sublattices B and W, see 
Fig.[2j The two spanning trees are dual, and from either one, the original dimer configuration can 
be entirely reconstructed. 




FIG. 2: (a) Dimer tiling; (b) spanning tree on the odd-odd sublattice or (c) even-even sublattice. 



re branches of the 



A directed path from a site x G B (x G W) to a site y G B (y G W) along t : 
black (white) spanning tree is the chemical path from x to y. As shown in 13t-ll7j|. the statistical 
properties of the chemical paths coincide with those of the loop erased random walk. Our main 
interest here is the return probability of the LERW: the probability that a chemical path starting 
at a given point on the spanning tree visits one of neighbouring sites of this point, say the right 
one, before going off to the root located at r. 

We consider now packed dimer tilings on the lattice £' obtained from the (2n — 1) x (2n — 1) 
lattice by removing three sites: the corner r, a grey site i and a black site j. The grey site i has 
two nearest neighbors G B, see Figj3l Kenyon made the observation (lemma 17 in 36(]) that 
if j was on the boundary of then the sites i\ and %2 would have to be in different components of 
the two-component black spanning tree, one component being rooted at r, the other at j. Indeed, 
i\ and 12 cannot be in the same component for otherwise the chemical path from i\ to 12 together 
with the bond of B linking i\ and 12 would form a loop enclosing an odd number of lattice sites of 
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(a) (b) (c) 

FIG. 3: (a) Dimer tiling with two monomers inside the lattice; (b) and (c) spanning trees corresponding to 
case (I). 

£ which cannot be fully covered by dimers. When j is not on the boundary, this result does not 
hold because one of the paths from %\ or from %2 may form a loop around j (so that one of i\,i2 is 
directed to r, the other to a loop winding around j), or else both i\ and 12 are directed to r, the 
chemical path between the two going around j. 

For our purposes, we will choose the site j as the left nearest neighbour of i\ on the sublattice 
B. Then a closed loop around j is impossible for the lack of space between j and %\ or i<i on the 
sublattice B. Three possibilities remain. 

(I) The first possibility is depicted in Fig. [3b: the path from i\ goes to j and the path from 12 
goes to r. By reversing the orientation of the first (red) path and inserting an arrow from i\ to 12, 
we obtain a long path from j to r passing through the bond (ii, ^2)- Depending on the bonds used 
by this path around i\, we have three possible configurations shown in Fig.[H 




FIG. 4: The tree possible path configurations for case (I). 

The first two are exactly those we had denoted by B and C in the previous section (up to a 
mirror transformation). The third one corresponds to all spanning trees on B which use the two 
bonds (J, i\) and (11,12). In terms of dimers on the original lattice C, it is the set of all dimer 
coverings with two forced dimers on j and i±. So the fraction of these is a local dimer-dimer 
correlation, computed by Fisher and Stephenson in Q], and equal to 1/47T in the limit of large 
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lattices. 

(II) The path from i\ goes to r and the path from 12 goes to j, see Fig. [5b. Like in case (I), we 
reverse the orientation the second path and insert an arrow form 12 to i±, to obtain a long path 
from j to r passing through the bond (i2,h)- It leads to the two possible local configurations A 
and B shown in Fig.[H 




A B 

FIG. 6: The two possible path configurations for the case (II). 

(Ill) The sites i\ and ii belong to the same component, rooted at r, and the chemical path 
between the two loops around j, almost completely enclosing the component rooted at j, see Fig. [7b 
and Fig. [8b- It implies that the dual graph on W contains one loop around the site j (and only 
one for the lack of space between i and j) . This loop can have the two orientations, as pictured 
in red in Fig. [7b and Fig. [8b. Rotating clockwise by ir/2 the red vertical arrow of the loop running 
in between i and j, we obtain a set of paths of type B or C, but now on the sublattice W (Fig.^J). 
The number of paths must be doubled since there are two equivalent orientations of the arrows 
along the loop. 

Collecting the various contributions from cases (I) to (III), we find how the dimer configurations 
on £' is related to the specific classes A, B, C of spanning trees on B. Dividing by the total number 
of dimer configurations on C, we obtain that the correlation P mm of two monomers at i and j is 
expressed in terms of the relative fractions of spanning trees of types A, B, C, in the limit of large 
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FIG. 9: Configurations of path for case (III). 



lattices, 



°mm = Mi + M 2 + M 3 = — + A + 4B + 3C. 

47T 



(3.1) 



Following [6(, the calculation of the correlation is easily carried out. One finds P mm = yielding 
a first relation for the three unknowns. 

In order to write a second relation, we repeat the previous calculation in which, in addition to the 
two monomers at i and j, we force a dimer in between them, like shown in Fig.llUb. Thus instead of 
a monomer-monomer correlation, we now consider a monomer-dimer- monomer correlation P m dm- 
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(a) 



(b) 



(c) 



FIG. 10: (a) two monomers; (b) two monomers and one dimer; (c) two dimers equivalent to case (b). 



All the steps above remain, except that case (I) becomes forbidden (the red loop in Fig. [3b would 
enclose an odd number of sites), and only one orientation of the loop in case (III) is allowed. 
Therefore we can write 



(3.2) 



Pmdm = M 2 + ^M 3 = A + 2B + C. 
The correlation P m dm is equivalent to have two fixed dimers, as shown in Fig.llOb. From |6], we 

find P m dm — g ■ 

Finally we notice the identity A = B as follows from the steps shown on Fig. [11] where we 
reverse the orientation of the loop, move the vertical up arrow to the horizontal right arrow and 
eventually apply a (diagonal) mirror transformation. 



A = 




f f f- 





= B 



FIG. 11: Three steps proving the equality of A and B. 



The two equations (|3.1j) and (|3.2p can then be solved, with the result 

32 4vr 

c-i— 5-. 

2vr 32 



(3.3) 
(3.4) 



Plugging these values back in the LERW return probability (j2.11j) yields P ret = 5/16, and in turn 
the ASM mean height (h) = 25/8. 



IV. PERSPECTIVE 



The present work raises (at least) two natural questions. We have shown how to compute 
the ASM mean height, a quantity so far thought to have a non-local interpretation in terms of 
spanning trees, in a purely local way, though specific local arrangements of dimers and monomers. 
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p c omp uting 2-site height correlations 



28 



It would be interesting to see whether this technique may he 
for heights larger or equal to 2 (since the others are known 

The second question is related to the LERW. Beyond the return probability, or passage probabil- 
ity to a nearest neighbour, one could ask for the passage probability to a second nearest neighbour 
(distance \/2). From numerical simulations, this latter probability appears to be close to the ra- 
tional value 2/9. Applying the present techniques to this case, and perhaps to the next few cases, 
would bring valuable results. 
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